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Abstract 

If the dark matter particle is a Majorana fermion, annihilations into two fermions 
and one gauge boson could have, for some choices of the parameters of the model, 
a non-negligible cross-section. Using a toy model of leptophilic dark matter, we 
calculate the constraints on the annihilation cross-section into two electrons and 
one weak gauge boson from the PAMELA measurements of the cosmic antiproton- 
to-proton flux ratio. Furthermore, we calculate the maximal astrophysical boost 
factor allowed in the Milky Way under the assumption that the leptophilic dark 
matter particle is the dominant component of dark matter in our Universe. These 
constraints constitute very conservative estimates on the boost factor for more 
realistic models where the dark matter particle also couples to quarks and weak 
gauge bosons, such as the lightest neutralino which we also analyze for some con- 
crete benchmark points. The limits on the astrophysical boost factors presented 
here could be used to evaluate the prospects to detect a gamma-ray signal from 
dark matter annihilations at currently operating lACTs as well as in the projected 
CTA. 
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1 Introduction 



The indirect search of dark matter has entered in the last few years into a new era 
of precision measurements. The Fermi Large Area Telescope (LAT) has measured with 
unprecedented accuracy the flux of cosmic electrons with energies between 7 GeV and 
1 TeV [1] and the flux of cosmic gamma-rays between approximately 20 MeV and 300 
GeV [2]. Besides, the satellite-borne experiment PAMELA has measured the cosmic 
antiproton flux and the antiproton-to-proton fraction between 60 MeV and 180 GeV [3], 
the electron flux between 1 and 625 GeV and the positron to electron fraction between 
L5 and 100 GeV [5]. 

The quality of the data collected by these experiments opens excellent opportunities 
to search for the annihilation or the decay of dark matter particles with masses in the 
electroweak range. Unfortunately, this search is hindered by the existence of large, and 
poorly known, astrophysical backgrounds which complicates in general the identification 
of a dark matter component in cosmic rays against the astrophysical backgrounds. 

This generic difficulty can be circumvented in some situations. The most notable 
one is when the dark matter particle annihilates or decays producing monoenergetic 
gamma-rays, which produces a very peculiar signal in the energy spectrum which cannot 
be mimicked by any known astrophysical source. Conversely, the non-observation of a 
gamma-ray line in the Fermi-LAT data sets very stringent constraints on the dark matter 
annihilation cross-section or decay rate into monoenergetic gamma rays P,[7j. 

On the other hand, when the dark matter is a Majorana particle, higher order pro- 
cesses such as XX ~^ f~^f~l can become important for some choices of the parameters of 
the model, despite being suppressed by the extra coupling and by the phase space fac- 
tor. In this case, gamma rays can be produced in the final state radiation, i.e. radiated 
off a nearly on-shell particle, or in the (virtual) internal bremsstrahlung, radiated 
off the internal bosonic line or from an off-shell f"^ particle [8], [9l [lO] . Whereas in the 
former process the energy spectrum of photons is rather featureless, and hence difficult 
to disentangle from the also featureless spectrum produced by astrophysical sources, the 
latter displays a very prominent bump which could be easily discriminated from the 
background. Several papers have recently analyzed the constraints on the annihilation 
cross-section into gamma-rays via internal bremsstrahlung from the Fermi-LAT mea- 
surements of the diffuse gamma-ray background JT] or the dwarf galaxy Segue 1 [12], 
from the H.E.S.S. observations of the globular clusters NGC 6388 and M 15 [13] and the 
dwarf galaxies Sculptor and Carina [11] , as well as from the MAGIC observations of the 
dwarf galaxy Willman 1 [15] . Furthermore, the prospects to observe this gamma-ray sig- 
nature at currently operating Imaging Air Cherenkov Telescopes (lACTs) was evaluated 
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in [ini [T7] and at the projected Cerenkov Telescope Array (CTA) in [161 [HI HH]- 

Due to the SU {2)lxU{1)y invariance of the Standard Model, the annihilation process 
XX ~^ f^f~1 is necessarily accompanied by the process XX ~^ f~ ^i which leads to 
the production of antiprotons through the fragmentation of the Z boson. Additional 
antiprotons can arise through the fragmentation of the fermions, provided they carry 
baryon number. Given the very stringent limits on the existence of an exotic component 
in the antiproton-to-proton fraction, it is important to analyze the constraints on the 
parameters of the model from the PAMELA measurements, in order to evaluate the 
prospects to observe the associated gamma-ray signal at present and future experiments. 

In this paper we will consider a toy dark matter model, where the dark matter particle 
is a Majorana fermion, singlet under the Standard Model gauge group, which couples 
to the electron lepton doublet and to a scalar doublet via a Yukawa coupling. Despite 
the dark matter particle being "leptophilic" the annihilation will necessarily produce 
antiprotons via electroweak bremsstrahlung processes. In Section [2] we will present the 
details of the model and we will calculate the cross-sections for the 2 — j- 3 processes 
involving two fermions in the final state, as well as a photon, a Z boson or a boson for 
different values of the parameters. In Section [3l we will calculate the limits on the cross- 
section for the processes involving weak gauge bosons from the PAMELA measurements 
on the antiproton-to-proton fraction. Finally, in Section H] we investigate implications 
for realistic models, in particular for several MSSM benchmark models that have been 
considered in the past in connection with internal bremsstrahlung. Various analytical 
expressions for 2 — 3 processes can be found in AppendixIXl 

2 Annihilations into two fermions and one gauge bo- 
son 

We consider a toy model, consisting in extending the particle content of the Standard 
Model by a Majorana fermion, x = X^-, which we assume to constitute the dark matter 
of the Universe, and a scalar field, rj. Their quantum numbers under the gauge group 
5t/(3), X SU{2)l X U{1)y are: 



Namely, x is a singlet under the gauge group while rj has the same gauge quantum 
numbers as the Standard Model Higgs. We further impose that the field rj carries electron 
lepton number Lg = — 1 in order to forbid its couplings to the quarks as well as to the 




(2.1) 
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second and third lepton generation^. With this particle content, the Lagrangian reads: 

r' p \ r \ r \ /'fcrmion , /'scalar /o 
L — i-SM + -L-x + '-r] + Mnt + Mnt • l-^-^J 

Here, £sm is the Standard Model Lagrangian which includes a potential for the Higgs 
doublet = m^^'f^ + iAi($^$)^. On the other hand and are the parts of the 
Lagrangian involving just the Majorana fermion x ^iid the scalar particle 77, respectively, 
and which are given by 

l^x = Ix^^^X - ^m^x'x , 

^ ^ ^ (2.3) 

where denotes the covariant derivative. Lastly, £fnt™^ ^"^^ ^I'i^^^^ denote the fermionic 
and scalar interaction terms of the new particles to the electron doublet and to the Higgs 
doublet: 

^r'"'^ = fx{L,^<y2V) + h.c. = fx{v,Lv' - CLV^) + h.c. , 

^scalar ^ _X,(^^ ^)(r^^^) - X,i^^^)i^^^) . 

We will assume that ml < and m2 > 0, therefore the minimization of the po- 
tential leads to (0°) = v^-m^/Ai, (r/°) = 0. The mass spectrum of the scalars can be 
straightforwardly calculated, the result being: 

m^o =ml + (A3 + X4)vlw ) (2-5) 
"^Ji = "^2 + ^svlw ■ 

Note that 77° is, like the charged component, a complex field that carries lepton number 
and contains the CP-even and CP-odd neutral components of the doublet 77. In order 
to guarantee the stability of the dark matter particle, we require that < m^o ,,±. For 
dark matter masses in the range 100 GeV- TeV, as we will assume here, this condition 
is automatically satisfied if m2 ^ Mz- 

The quartic couplings in the scalar potential are subject to constraints from demand- 
ing the absence of unbounded from below directions in the potential. These conditions 
were derived in [20] : 

Ai>0, A2>0, A3 > -(A1A2)'/' , 

(2.6) 

A3 + A4 > -(AiAs)^/^ . 



similar model, with lepton number softly broken, has been considered in [19] in connection to 
neutrino masses. 
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These four conditions can be simultaneously fulfilled for both signs of the coupling A4. 
Therefore both orderings < m^o and > m^o are theoretically viable, which as 
we will see has implications in determining which is the dominant annihilation channel 
of the dark matter particle. 

The Majorana dark matter particle x can be produced thermally in the early uni- 
verse, and the correct abundance could be achieved, for a given dark matter mass, by 
adjusting the coupling constant / [21] (neglecting coannihilations, and assuming a freeze- 
out temperature T/o = m£)A//20 and g^^iTjo) = 106.75, see e.g. [22]): 



^^o.iif^Vr^^v 



QoAih' ~ 0.11 



/ J VlOOGeV/ 
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(2.7) 

In principle, this would fix also the absolute strength of any signals relevant for dark 
matter detection. However, since we are here mostly interested in a comparative analysis 
of different indirect dark matter signals, and since the production could be influenced by 
additional heavier degrees of freedom not included in the toy model or by non-thermal 
effects, we will not a priori restrict the value of / in the following (unless stated otherwise) 
but rather assume that the particle x constitutes the dark matter in the present universe 
with the correct abundance. 

The observational prospects for detecting energetic particles produced in annihi- 
lations of dark matter particles are determined by the velocity-weighted annihilation 
cross-section. It can be expanded into s- and p-wave contributions, 

av = a + bv^ + 0{v^) . (2.8) 

For typical relative dark matter velocities v < 10~'^c within the Milky Way halo [23], the 
p-wave contribution bv"^ is usually strongly suppressed. 

The lowest-order annihilation channels are the two-to-two processes XX ~^ and 
XX ~^ which proceed by exchanging an ?7^-particle or an ?7*^-particle, respectively. 
Because of the Majorana nature of x? these processes can proceed via the t- and u- 
channel. The interference of these diagrams leads to the well-known helicity-suppression 
of the s-wave contribution to the two-to-two annihilation cross-section, a oc rnl/m\)j^,p 
In the limit mg = 0, the cross-section for the channel xX ~^ is given by [21] 

/4 l+m^±/m^ 



For XX ~^ one has to replace m^± H- m^o. Since the s-wave contribution is helicity 
suppressed, and the p-wave contribution is suppressed by v"^, the annihilations xX ~^ 
and XX ~^ are not efficient in the Milky Way today. 



Additional annihilation channels arise at higher orders in perturbation theory, e.g. 
the two-to-three process XX ^ 1^^- Due to additional couplings and phase-space, one 
would naively expect a suppression factor of order aem/T^ ~ 2 ■ 10"'^ compared to the 
two-to-two process. As is well-known, the emission of a soft or collinear photon off the 
final electron (or positron) is logarithmically enhanced ~ aem/7r[ln(mg/(4m|,j^^))]^, and 
leads to a model-independent spectrum characteristic for final state radiation (e.g. [23]). 

More importantly, the emission of a photon can also lift the helicity suppression. The 
reason is that the three-body final state 7ee can simultaneously accommodate a left- 
handed electron, a right-handed positron, and possess zero angular momentum. There- 
fore, the s-wave contribution is not helicity suppressed. Consequently, the annihilation 
channel xX ~^ 7^6 can even dominate over the lowest-order one, because the absence of 
helicity suppression more than compensates the suppression factor aem/T^ [H]- 

Note that this feature is an important difference to the soft or collinearly emit- 
ted bremsstrahlung photons, which yield only a comparably small contribution. The 
copious emission of hard and non-collinear photons has been called 'virtual internal 
bremsstrahlung', because the emission can be viewed as a part of the hard process. We 
will refer to this process simply as internal bremsstrahlung (IB) in the following. Typi- 
cally, the spectrum of the produced photons is peaked at high energies, of the order of 
the dark matter mass. Therefore, these photons lead to a characteristic gamma-ray signal 
that can be potentially observed by present and future gamma-ray telescopes p ^ fTT l fTS]. 

In addition to the two-to-three process xX ~^ 7^6 producing a photon, SU{2)l x 
f/(l)y-invariance inevitably leads to the presence of the two-to-three process xX ~^ Zee, 
provided the channel is kinematically allowed, i.e. m^i^ > Mzjl- As for the pho- 
ton, one can discriminate between ordinary final state radiation and (virtual) internal 
bremsstrahlung. In the former case, the Z-boson is radiated off an on-shell electron. Al- 
though this process cannot lift helicity suppression, it is subject to a logarithmic enhance- 
ment and can be important for heavy dark matter particles, rtiDu ^ TeV [25 | [26 | [27 t [28] . 

On the other hand, internal electroweak bremsstrahlung describes the emission of 
a Z-boson off the internal line, or from an off-shell electron. This process can lift the 
helicity suppression, like the photon. It is the main purpose of this work to study these 
electroweak IB processes, for which one can expect a similar enhancement as for the 
electromagnetic IB process xX ~^ 1^^- The further decay and fragmentation of the Z- 
boson produces in particular protons and antiprotons, in equal amounts. In addition, 
the decay products lead to additional photon, electron and positron as well as neutrino 
production. 

In view of the stringent observational constraints on the cosmic flux of antiprotons, 
it is important to study the process xX ~^ Zee, in particular when considering model 
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Figure 1: Spectrum of photons as well as Z- and W^-bosons produced by internal 
bremsstrahlung in dark matter annihilations for mrtu = 300GeV and m^± = m^o = 
330GeV. 

parameters that lead to a strong gamma-ray signal from electromagnetic IB. Within the 
toy model considered here, also the electroweak processes XX Zuu, xx ~^ Weu, and 
XX ~^ Weu contribute to the production of antiprotons. 

The energy distribution of H^-bosons produced by IB in dark matter annihilations 
can be represented by the dimensionless differential cross-section 

dNw _ 1 / vdaixx^Weu) vda{xX^Weu) \ 

dlnE av{xx^ee) \ d\nE dlnE ) ' ^ ' ' 

which is defined analogously to the corresponding quantity dN^/dlnE = av{xx ~^ 
ee)~^ vda{xx ~^ '-fee) /dlnE for electromagnetic IB. The normalization to the two-to- 
two cross-section (tv{xX ~^ ^c) is chosen as in [8] for electromagnetic IB for convenience. 
Similarly, one can also define dNz/dlnE. 

Our analytic results for the s-wave contribution to the differential cross-section can 
be found in appendix |X1 They are applicable for arbitrary mass spectrum, except for 
neglecting the electron mass mg = 0, in the relevant range Miy/2 < rriDM < m^±,m^o, 
in particular for masses m^^ nearly degenerate with tudm and also for m^± ^ m^o. The 



numerical results shown in the following include also velocity-suppressed contributionqj. 



■^We have used CALCHEP [3SJ [35] for parts of the analytical as well as for numerical computations. 
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The W, Z, and photon spectra from IB are shown in Fig.[Tl The massive gauge 
bosons are produced with a hard spectrum that is peaked near the maximal energy, 
similar to electromagnetic IB. However, the maximal energy is slightly higher, Emax = 
mj:,M + M'^^y^-^/^AmcM)- Another difference is the lower cut-off at the gauge boson mass. 
Since the shape of the spectra near the peak are nevertheless very similar, we will use the 
total cross-sections in order to illustrate the relative importance of the various processes 
in the following (we use the full spectra in our numerical analysis, however). 

The relative strength of electromagnetic IB, XX ~^ 7ce, and the electroweak IB 
process xX ~^ Zee, for dark matter masses far above the threshold rriDM ^ Mzjl-, 
is simply related by the relative couplings of (left-handed) electrons to photons and to 
Z-bosons, respectively, 

ai;(xx^^ee) : (Tt;(xx ^ 7ee) = cot2(2^,4/) = 0.41 for m^M > . (2.11) 

This is because both processes are mediated by the same particle, namely the charged 
component rf" of the doublet. For the other electroweak IB processes, also the neutral 
component rf of the doublet appears in the relevant Feynman diagrams. Therefore, their 
branching ratios depend on the mass m^o, or, more specifically, on the mass splitting 



The case m^o = m^± 



Let us first assume that the components of the doublet are degenerate, m^o = m^±, 
i.e. A4 = in the scalar potential. Since both mediating particles are degenerate, all 
branching ratios are simply determined by the relative coupling constants for m£,M ^ 
Mz/2, 



(^v{xX Zuv) : av{xx ^ 7ee) = —^^^^ = 1-41 ^ M. 



1 _ ^ for TTiDM > '-7^-, = m^± 

(2.12) 



(^v{xx Weu) : av{xx ^ 7ee) = —^Jq^) = 4.32 2 



Here crvixx ~^ Weu) = crvixx ^ ^ ez/) + avixx ~^ W^ev). 

The full result for the ratio of the two-to-three cross-sections to the two-to-two cross- 
section (jv{xx ~^ se) is shown in Fig.[2] for m^o = m^±, as a function of the dark 
matter mass. Clearly, the two-to-three processes dominate over the two-to-two annihi- 
lation channel, i.e. all ratios are much larger than one. Furthermore, one can observe 
that the branching ratios obey the simple relations given above for dark matter masses 
"^DA/ ^ 200GeV. For dark matter masses in the range 45GeV < ^ 200GeV, the 



^This case has also been recently studied in Refs. [43l |44] 
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Figure 2: Ratio of three-body and two-body annihilation cross-sections, for electro- 
magnetic IB, (jv{xx ~^ 7ee)/crf(xx ~^ ^e), and for the electroweak IB channels 
XX ~^ Zee, XX ~^ Zvv and XX ~^ Wev. The latter denotes the sum of W~ev and W^ev. 
The ratio of the scalar doublet mass to the dark matter mass is fixed to m^± /m^A/ = 1-1 
(left) and 1.5 (right), and here we assume that both components of the doublet have 
identical mass, m^o = m^±. For the relative dark matter velocity we use v = 10~^c. 

branching ratios for electroweak IB become kinematically suppressed, but are still size- 
able for rriDM ^ lOOGeV. Since the energies we are interested in, i.e. at which lACTs 
are sensitive to observe a gamma-ray signal, typically lie above ~ lOOGeV, we generi- 
cally expect an 0{1) branching fraction into electroweak gauge bosons in addition to the 
gamma signal from IB. 

In Fig. (31 the ratio of two-to-three and two-to-two cross-sections are shown as a 
function of the mass m^± of the particles rj^ and 77*^ that are exchanged in the annihilation 
process, assuming m^± = m^o. Similar to electromagnetic IB, the electroweak IB is 
strongly enhanced when m^± is nearly degenerate with rriuM- The IB processes dominate 
over the two-to-two channel as long as both masses are of the same order of magnitude, 
1 < mrj± /moM < 5. In the limit ■mrj±/m£,M ^ 1, we find that the s-wave electroweak 
IB cross-sections scale like av{2 — )■ 'i)\s~wave oc {mr)M/mrj±)^. This is the same scaling 
compared to electromagnetic IB. Since the two-to-two cross-section scales like av{2 — )■ 
2) oc {m£)M/m^±)^, all the IB processes become suppressed when the mediating particles 
are very heavy, roughly m^± /rrij^M ^ 5. In that case the conventional final state radiation 
of soft and collinear photons would be the dominant production mechanism, which scales 
also with {rriDM /'mr)±Y 1 but cannot lift the helicity suppression. However, there is a large 
and generic region of parameter space where electroweak and electromagnetic IB are the 
dominating annihilation channels. 
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m„± /niDM 



Figure 3: Ratio of three-body and two-body annihilation cross-sections for electromag- 
netic and electroweak IB, as a function of the doublet mass m^±. The dark matter mass 
is rriDM = 300GeV, and m^o = m^±. Similarly to electromagnetic IB, the strongest 
enhancement of electroweak IB occurs when the dark matter mass rriDM and the mass 
m^± of the mediating particle are nearly degenerate. The IB annihilation channels domi- 
nate for m^± < 5m£,M, i-e. when the masses are of comparable order of magnitude. For a 
strongly hierarchical spectrum, m^± > 5mDM, the two-to-two channel as well as collinear 
electromagnetic final state radiation dominate. For the relative dark matter velocity we 
use V = 10~^c. 
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The case m^o ^ m^± 

Due to electroweak symmetry breaking, one expects in general that the neutral and 
charged components of the SU{2)l doublet t] acquire different masses. The mass dif- 
ference is controlled by the order parameter vew = 174GeV of electroweak symmetry 
breaking. In particular, within the toy model one finds that m^o — fn'^± = X4v'^^r, where 
A4 is a coupling in the scalar potential that is generally non-zero. Also, for example 
within the CMSSM, the components of the slepton doublet are non-degenerate at the 
electroweak scale. Therefore, it is important to discuss this case. For the toy model con- 
sidered here, this case has not been discussed in the literature in connection with IB to 
our knowledge. 

Clearly, the cross-sections for xx ~^ 7^^^ XX ~^ are not affected by a mass- 
splitting, since they only involve the rj'^ particle. Furthermore, the channel xX ~^ '^'^^ 
just involves the rf particle, and therefore (neglecting lepton masses) 

/ _|_ \ 

av{xx Zuu) = ^^^2^) (TvixX ^ Zee)\„^^^^m^, ^ 1.41 av{xx ^ 7ee) , 

(2.13) 

where the latter estimate is valid for m^i ^ Tn^i^j ^ Mz/2. This is precisely the be- 
haviour expected from simple scaling arguments. In particular, it results in an additional 
enhancement of electroweak IB when m^o < m^± . Since we are primarily interested in the 
case of a potentially strong gamma signal, we will assume that m^o > m^±, i.e. A4 > 0, 
in the following. 

The most interesting process is the channel xX ~^ Wev. The first reason is that this 
is the electroweak IB process which yields the largest contribution within the toy model, 
and is thus the crucial one for antiproton production. The second reason is related to 
the fact that both the neutral and charged components iff and rj"^ simultaneously appear 
in the corresponding Feynman diagrams. In particular, the vertex rf'q'^W^ features a 
coupling to both the transversally polarized as well as longitudinally polarized VF-bosons, 
Wt and Wl. The latter is related to the coupling to the Goldstone bosons that give 
mass to the H^-bosons, according to the equivalence theorem. This coupling is given by 

£D-4 "^""^^% Vg-+h-c.. (2.14) 

Mw ^ 

Thus, the coupling is non-zero only in the presence of a mass splitting. This means 
that in the case m^o = m^± discussed above only transversely polarized 1^-bosons are 
produced, xX ~^ Wt^v. When taking a mass-splitting into account, a new channel opens 
up, namely xX ~^ Wlcv. This feature represents a genuine difference compared to the 
IB of photons. Note that the upper vertex actually contributes in Feynman gauge, while 
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Figure 4: Ratio of three-body and two-body annihilation cross-sections (7v{xx ~^ 
Weu) /av{xx ~^ ce) for the electroweak IB channels into W-bosons {Weu = W~ev + 
W^ep). The cross-section avixx ~^ Weu) depends sensitively on the mass-splitting 
of the components of the doublet. For the previously considered case where 
: m^± (dot-dashed line), only transversely polarized W-bosons can be produced. 
For comparison, we show the cross-section when taking a mass-splitting into account 
(solid and dashed lines). In that case also longitudinally polarized W-bosons can be pro- 
duced. In particular, the choice m^o — '^^i = ^iiy with vew = 174GeV characterizes 
a typical splitting expected from electroweak symmetry breaking. The other parameters 
are chosen as in fig. [21 



it is implicitly taken into account in unitary gauge. We have checked that our results 
obtained in both gauges agree with each other. 

The annihilation XX ~^ WlCV is important for two reasons: first, we find that for 
a generic mass splitting m^o — 'm1^± = '^iv\^r with A4 ~ the annihilation into 

longitudinal M^-bosons significantly enhances the total cross-section. This can be seen 
in Fig.m where the case with mass splitting is compared to the degenerate case. The 
enhancement is important for dark matter masses up to m/jjv/ ~ 1 ~ 2TeV. A simple 
analytic estimate is possible for m^, ^ ttidm ^ 



1 



sm 



4.32 



5 (mjo - 



2m2, 



m^o -|- . 



1 + A 



/aooGevy 

V f^DM J 



(^v{xx 7ee) 
(^v{xx -> 7ee) ,(2.15) 



where the second term in the square bracket corresponds to longitudinal Vl^-bosons. More 
accurate expressions can be found in appendixRl 

Second, also the spectrum of the produced longitudinal Wl differs from the spectrum 
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Figure 5: Spectrum of ly-bosons produced by IB in dark matter annihilations for rriDM = 
300GeV, m^± = 330GeV and different choices of m^o. For m^o ^ 7Ti^±, also longitudinal 
VT-bosons are produced, which can have the maximal energy Emax = itt-dm + / AmuM 
and therefore lead to a harder spectrum. 
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Figure 6: Ratio of three-body and two-body annihilation cross-sections (Jv{xx ~^ 
Wev)/ av{xx ^ ^e) for the electroweak IB channel into W-bosons, for different choices 
of the mass-splitting m^o — "i^± , as a function of the charged scalar mass m^± . The other 
parameters are chosen as in Fig.|3l 



of the transverse Wt considered before. The reason is that for xx ~^ Wieiy, it is possible 
that the e and u are emitted in one direction, while the Wl is emitted in the opposite 
direction. This configuration supplies the Wl with the maximal possible energy Emax = 
ruDM + / ArriBM- For a transversal Wt, this configuration is forbidden by angular 
momentum conservation (up to a tiny helicity-suppressed contribution). Therefore, the 
spectrum of VF-bosons obtained when taking a mass splitting into account is significantly 
harder than the spectrum obtained in the case m^o = m^± . This behaviour can be clearly 
seen in Fig.[5l and is also manifest in the full analytic result presented in appendixRl 

In Fig.[6l the scaling of the cross-section (Jv{xX ~^ Wev) with m^± is shown, when 
keeping the mass-squared difference m^o — = ^^'^ew dark matter mass 

fixed. As before, the s-wave contribution scales with {m£)M/''^r]±)^- We find that the 
enhancement of the cross-section due to emission of longitudinal Wl persists for large 
m^± ^ ttldm, even when keeping m^o — ^'^± constant, as appropriate. 
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3 Antiproton flux and observational constraints 



We consider a distribution of dark matter particles in the Milky Way given by p{'r), 
where r denotes the position of the dark matter particle with respect to the center of 
our Galaxy. In order to estimate the uncertainty in the predictions stemming from our 
ignorance of the actual dark matter distribution in our Galaxy, we will calculate the 
predictions for three different halo profiles. Concretely we will assume for simplicity a 
spherically symmetric distribution with a radial dependence given by the Isothermal 
profile 

p{r) = , , f', S2 > (3-1) 
1 + [r/rsY 

the Navarro- Frenk- White (NFW) profile [23 ED] 



and the Einasto profile [311 132] 



p(r) = exp <^ 



r_2 




1 


[ a 







(3.3) 



In these expressions, the scale radius is = 4.38, 24.42 and 28.44 kpc, respectively, 
while a = 0.17 for the Einasto profile. Besides, the parameter ps is adjusted in order to 
yield a local dark matter density p{rQ) = 0.39 GeV/cm^ [33] with = 8.5 kpc being the 
distance of the Sun to the Galactic center and is given, respectively, by ps = 1.86, 0.25 
and 0.044 GeV/cml 

Two dark matter particles at the position r can annihilate producing antiprotons at 
a rate per unit of kinetic energy and volume given by: 

where (av) f is the thermally averaged cross-section multiplied by the velocity in the 
annihilation channel / and dNp /dT is the energy spectrum of antiprotons produced in 
that channel per unit of kinetic energy. 

The weighted antiproton spectrum {av) f dNp / d\n T / {av)^^^ee, normalized to the 
two-to-two cross-section is shown in Fig. [3 The left part shows the spectrum resulting 
from the various electroweak IB bremsstrahlung processes and the right part the spec- 
trum resulting from VT-production for various values of the scalar mass-splitting. We 
have obtained the spectrum using the event generator PYTHIA 8.1 [31] interfaced with 
CALCHEP [SHIES]. 
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Figure 7: Spectrum of antiprotons (before propagation) resulting from the decay and 
fragmentation of the W- and Z-bosons produced in electroweak IB processes. Parameters 
are itidm = 300GeV, m.;,± = 330GeV. The left plot shows the individual contributions 
for m^o = m^± , and the right one the W contribution for different choices of the mass- 
splitting mlo - mK . 



After being produced at the position r, antiprotons propagate through the Milky 
Way in a complicated way before reaching the Earth. Following [H7], we will describe an- 
tiproton propagation by means of a stationary two-zone diffusion model with cylindrical 
boundary conditions. Under this approximation, the number density of antiprotons per 
unit kinetic energy, fp{T,f,t), approximately satisfies the following transport equation: 

df- 

= ^ = V ■ {K{T, r}Vfp) - V ■ (K(r)/p) - 2h6{z)T,^Jp + Q{T, r) . (3.5) 

The boundary conditions require the solution fp{T, r, t) to vanish at the boundary of the 
diffusion zone, which is approximated by a cylinder with half-height L = 1 — 15 kpc and 
radius i? = 20 kpc. 

The first term on the right-hand side of the transport equation is the diffusion term, 
which accounts for the propagation through the tangled Galactic magnetic field. The 
diffusion coefficient K{T, r) is assumed to be constant throughout the diffusion zone and 
is parametrized by: 

K{T) = Kof3n' , (3.6) 

where P = v/c and TZ is the rigidity of the particle, which is defined as the momentum 
in GeV per unit charge, TZ = p{GeV)/Z. The normalization Kq and the spectral index 
6 of the diffusion coefficient are related to the properties of the interstellar medium and 
can be determined from the flux measurements of other cosmic ray species, mainly from 
the Boron to Carbon (B/C) ratio [38] • The second term is the convection term, which 
accounts for the drift of charged particles away from the disk induced by the Milky 
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Model 


6 


Ko (kpcVMyr) 


L (kpc) 


K (km/s) 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX 


0.46 


0.0765 


15 


5 



Table 1: Astrophysical parameters compatible with the B/C ratio that yield the minimal 
(MIN), median (MED) and maximal (MAX) flux of antiprotons. 



Way's Galactic wind. It has axial direction and is also assumed to be constant inside the 
diffusion region: K;(r) = Vc sign(z) k. The third term accounts for antimatter annihilation 
with rate Fann, when it interacts with ordinary matter in the Galactic disk, which is 
assumed to be an infinitely thin disk with half-width h = 100 pc. The annihilation rate, 
is given by 

Tann = (^H + 4:^/^nne)a;;'^Vp . (3.7) 

In this expression it has been assumed that the annihilation cross-section between an 
antiproton and a helium nucleus is related to the annihilation cross-section between an 
antiproton and a proton by the simple geometrical factor 4^/'^. On the other hand, tt-h ~ 
1 cm~'^ is the number density of Hydrogen nuclei in the Milky Way disk, riHe ~ 0.07 riH 
the number density of Helium nuclei and app^ is the annihilation cross-section, which is 
parametrized by [39] : 

ann/^^N _ / 661 (1 + 0.0115 T-^'^^^ - 0.948 r°-0151) ^^g^^^ ^ rj^ ^ 15 5 Q^y ^ 

Lastly, Q{T,f) is the source term of antiprotons from dark matter annihilations, which 
was defined in Eq. (13.41) . In this equation, energy losses, reacceleration effects and non- 
annihilating interactions of antimatter in the Galactic disk have been neglected. The 
ranges of the astrophysical parameters that are consistent with the B/C ratio and that 
produce the maximal, median and minimal antiproton fluxes are listed in Table [H 

Finally, the interstellar flux of primary antiprotons at the Solar System from dark 
matter annihilation is given by: 

^f{T) = ^fp{T,r^), (3.9) 



where v is the antiproton velocitjo. 

However, this is not the antiproton flux measured by experiments, which is affected 
by solar modulation. In the force fleld approximation [H] the effect of solar modulation 



*We have cross-checked our results using the interpolating functions presented in 
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Figure 8: Antiproton-to-proton ratio obtained from electroweak IB, for dark matter 
masses lOOGeV (upper left), 300GeV (upper right), ITeV (lower left) and 3TeV (lower 
right). Also shown are the PAMELA data for the p/p ratio [3]. The two lines corre- 
spond to the signal allowed at 95%C.L. (dashed), and to signal+background (sohd), 
respectively. The antiproton background is taken from |47j . 

can be included by applying the following simple formula that relates the antiproton flux 
at the top of the Earth's atmosphere and the interstellar antiproton flux |42j : 

where Tjs = Ttoa + 't'F^ with Tjs and Ttoa being the antiproton kinetic energies at the 
heliospheric boundary and at the top of the Earth's atmosphere, respectively, and (f)p 
being the solar modulation parameter, which varies between 500 MV and 1.3 GV over the 
eleven-year solar cycle. Since experiments are usually undertaken near solar minimum 
activity, we will choose (pp = 500 MV for our numerical analysis. 

In order to obtain constraints on the dark matter annihilation cross-section, we com- 



18 



pute the antiproton-to-proton ratio 

n + $bkg 

' 1 ' ^ (3.11) 

using the proton flux of [16]. Furthermore, in order to obtain a conservative exclusion 
bound we adopt the minimal value for the antiproton background as discussed in |47j . 

In Fig.|H]we show the resulting p/p ratio for various dark matter masses, together 
with the PAMELA data [3]. The data can be described consistently by the background 
only. Therefore, the measurements can be used to set an upper limit on the contribution 
to the antiproton flux originating from dark matter annihilations, which in turn 
yields an upper bound on the corresponding annihilation cross-section. 

We show in Fig.[9]our results for the upper limits on the cross-sections of the indi- 
vidual electroweak IB processes obtained from the PAMELA p/p data [3] using a x^-test 
at 95%C.L. In order to take the astrophysical uncertainties into account, we compute 
limits for the three propagation models and three halo profiles discussed above. How- 
ever, we emphasize that the limits on the individual channels are rather robust from 
the particle physics side. The reason is that they just depend on the shape of the an- 
tiproton spectrum, which turns out not to depend strongly on e.g. the scalar masses, as 
long as m^i/niDM ~ For the MED propagation model the limits lie in the range 

av < 10~^^cm^/ sec for itidm = lOOGeV and 10"^'^cm'^/ sec for uidm = ITeV for all the 
electroweak IB processes. 

In order to compare our results with other constraints, like e.g. the relic density, it 
is convenient to convert the exclusion limits on the cross-sections into constraints on 
the parameter space of the toy model. All the cross-sections relevant for dark matter 
annihilation are proportional to f^, where / characterizes the coupling strength of the 
dark matter particle to the Standard Model. Thus, the exclusion limits can be converted 
into upper limits on /, by using the total cross-sections derived in Section[21 

The resulting upper limits on the coupling / depend strongly on the ratio m^± /m^iM, 
because of the enhancement of the total IB cross-sections when mri± /ttidm 1. The 
results for several values of the ratio are shown in Fig.[TDl For comparison, we also show 
the coupling for which thermal freeze-out would yield the correct relic density within the 
toy model [21], as well as the exclusion bounds scaled with hypothetical boost factors 
BF in the Milky Way. For BF = 1, the upper bounds on the coupling from antiproton 
constraints are compatible with the couplings that yield the correct relic density, for all 
dark matter masses and all values of m^± /rriDM considered in Fig. [101 For sub-TeV dark 
matter masses, the antiproton constraints become relevant for very large values BF ~ 10^ 
of the boost factor, when assuming that x '^as produced thermally and constitutes the 
dominant component of dark matter. In that case, the antiproton constraints arising 
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Figure 9: Upper bounds on the electroweak IB cross-sections av{2 — )■ 3) obtained from 
the PAMELA data [3] of the cosmic antiproton-to-proton ratio. The plots correspond to 
the constraints on the individual electroweak IB processes XX ~^ Wev (where Wev = 
W'eu + W'^eu, upper), and xX ~^ Zee (lower) at 95%C.L. The bounds on xX ~^ Zvv 
are identical because the spectrum of Z bosons per annihilation is the same as for Zee. 
We have used m^±lraY)M = 1-5 and m^o — rn'^± = v%^r. Dashed lines correspond to 
the MIN, solid to MED and dotted to MAX propagation models. The bounds obtained 
assuming an isothermal dark matter profile are shown in red, NFW in black, and Einasto 
in blue. 
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10^ 10-^ 10' 10^ 10'' 10'' 

niDM [GeV] niDM [GeV] 

Figure 10: Upper bounds on the dark matter coupling / inferred from the antiproton 
constraints on the electroweak IB processes (95%C.L.) for m^o ""^^i = '^'ew ^^"^ different 
choices for the ratio m^± jnij^M = 1-01, 1.1, 1.5, 3.0. Also shown is the coupling for which 
thermal freeze-out would yield a relic density compatible with WMAP within the toy 
model [21] (dot-dashed), as well as the effect of a hypothetical boost factor BF in the 
Milky Way on the bound. The bounds correspond to the MED propagation model. 
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Figure 11: Upper bounds on the boost factor BF in the Milky Way inferred from the 
antiproton constraints on the electroweak IB processes (95%C.L.) for m^,, — 'ni'^± = v%y^ 
and different choices for the ratio m^± /rriDM = 1.01,1.1,1.5,3.0. The couphng / is 
fixed for each dark matter mass from requiring a rehc density from thermal production 
nh'^ = 0.11, see Eq.(^. 

from electroweak IB may be considered as conservative upper bounds on the size of the 
coupling /. 

Nevertheless, the constraints can be used to evaluate the possibility of observing 
signals from dark matter annihilations in other channels, in particular gamma rays or 
positrons, with present or planned instruments. For example, it has been argued that 
the positron excess observed in the PAMELA data [U [5], could be explained by dark 
matter annihilations into leptons. In particular, an enhancement of the annihilation into 
leptons due to internal bremsstrahlung has been considered [IH]. Note that in this context 
boost factors of the order ~ 10^ have been invoked in the literature. Despite a lack of 
explanation for astrophysical boost factors of this size, we note that the constraints 
coming from antiproton production due to electroweak IB independently restrict this 
interpretation of the PAMELA positron excess. Also, since the production mechanism 
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™dm [GeV| moM [GeV] 

Figure 12: Upper bounds on the boost factor BF in the Milky Way inferred from the 
antiproton constraints on electroweak IB processes (95%C.L.) for mfj±/m£)M = 1-1 (left) 
and 1.5 (right) and different choices for the mass-splitting m^o — ""^^i- The coupling / is 
fixed for each dark matter mass from requiring a relic density from thermal production 
Qh^ = 0.11, see Eq. (12. 7p . The bounds correspond to the MED propagation model and 
the NEW halo profile. 

of dark matter is unknown, it is worthwhile to investigate the viability of observably 
large dark matter annihilation cross-sections. In this sense, the boost factor can serve as 
a reference value in order to compare the cross-relations between prospects in different 
annihilations channels. 

Detecting a gamma-ray signal from internal bremsstrahlung at MAGIC II and CTA, 
for instance, typically requires boost factors in the target (e.g. the dwarf galaxies Draco 
or Willman 1) of the order of BF > 10^ — 10^ according to Ref. [16]. We stress that 
antiprotons from electroweak IB indeed yield relevant constraints for boost factors in the 
Milky Way of this size. 

If one assumes that the coupling / is fixed such that the correct relic density is 
produced in the early universe, then it is possible to obtain upper hmits on the boost 
factors BF that are compatible with antiproton constraints. The corresponding exclusion 
plots are shown in Eig.[TTl again for several values of the ratio m^±/m£)M, and assuming 
a generic scalar mass-splitting m^o — fn'^± = v^w- The results obtained when assuming 
m^o = m^± and m^o — m^± = 2f|;^, respectively, are shown in Eig.[T2]for comparison. 
We find that, for the case mrj±/mDM ~ 1, which yields a particularly strong gamma (and 
possibly positron) signal from electromagnetic IB, also the constraints from electroweak 
IB become quite stringent in comparison, of the order BF < 10^ — 10'^. 
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Model 






BF (p/p) 






[GeV] 


MIN 


MED 


MAX 






2 ^ 2/3(2 ^ 3) 


2 ^ 2/3(2 ^ 3) 


2 ^ 2/3(2 ^ 3) 


BM2 


453 


< 4.4-10^ (9.6 ■ 10^) 


< 5900 (1.3- 10^) 


< 2700 (5.8 • 10^) 


BM3 


234 


< 1.0-10^ (9.9-10^) 


< 1500 (1.3-104) 


< 660 (6100) 


BMJ' 


316 


< 2400 (2.6 ■ 10^) 


< 330 (3.5 ■ 10^) 


< 150 (1.6-104) 


Bur 


141 


< 82 (5.1 • 10^) 


< 11 (6900) 


< 5.0 (3100) 



Table 2: Upper limits on the boost factor BF in the Milky Way obtained from the 
PAMELA p/p data [3] for several MSSM benchmark points (see [U [16] for details) at 
95%C.L. We have taken into account the antiproton production from 2 — 2 processes 
with hadronic final states as well as 2 — t- 3 processes that are enhanced due to electroweak 
IB. The values in parentheses are the upper limits on the boost factor that would be 
obtained when taking only the 2 — 3 processes into account. The bounds correspond to 
the MIN, MED and MAX propagation models and the NEW halo profile. 

4 Implications for realistic models: the case of the 
lightest neutralino 

In the previous section we have derived bounds on the 2 — )■ 3 annihilation cross-section 
and the corresponding boost factor for a toy model of leptophilic dark matter which only 
couples to the electron doublet and an extra scalar doublet. However, in more realistic 
models the WIMP dark matter particle also couples to the quarks and the gauge bosons. 
As a consequence, the antiproton production in 2 — )■ 2 processes can be significant or 
even dominate over the production in the fragmentation of the weak gauge bosons from 
the 2 — )■ 3 processes which we have analyzed in this paper. It is then apparent that the 
bounds derived in this paper for the boost factors can be used as conservative bounds 
for any realistic model. 

In order to estimate how are these conservative bounds modified when consider- 
ing realistic models, we have studied the very popular scenario where the dark matter 
particle is the lightest neutralino. We have analyzed the antiproton production in the 
supersymmetric benchmark points BM2, BM3, BMJ', BMP (using [l9l[50]); the points 
BM2 and BM3 lie in the coannihilation region, BMJ' in coannihilation tail and BMP 
in the bulk region. All of them are characterized by a large cross-section for the inter- 
nal electromagnetic bremsstrahlung processes. Eurthermore, BM3, BMJ' and BMP were 
used in [16] to evaluate the prospects to observe a gamma-ray signal in dark matter 
annihilations at MAGIC II or at the projected CTA. Our results are shown in Table [21 
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where we list the upper hmit on the boost factor in the Milky Way from requiring no 
excess in the PAMELA p/p data at the 95% C.L., assuming the NFW halo profile and 
for the MIN, MED and MAX propagation models. We also show in parentheses the upper 
bound corresponding to just considering the antiproton production in the electroweak 
bremsstrahlung. From the table it follows that in realistic models the upper bounds 
on the boost factor are at least one order of magnitude stronger than the conservative 
bounds derived in this paper, and in some cases even three orders of magnitude stronger. 

It is interesting to compare our upper bounds on the boost factor in the Milky way 
with the lower bounds on the boost factor derived in [16] for the Draco dwarf spheroidal 
galaxy, from requiring the detection of a gamma-ray signal with a 5(T significance at 
MAGIC II or at the projected CTA with 30 hours of observation. The observation of 
such a signal at MAGIC II requires a boost factor in Draco larger than ~ 10^, and at 
the CTA larger than ~ 10^. In contrast, we obtain that the non-observation of an excess 
in the PAMELA p/p data implies, for the same benchmark point, a boost factor in the 
Milky Way smaller than ~ 10^ - 10^ for the MED propagation model. Therefore, the 
observation of a gamma-ray signal from Draco at MAGIC II will imply the existence of a 
boost factor in Draco much different to the boost factor in the Milky way. Furthermore, 
the supersymmetric benchmark point BM3 has been discussed in [18] as a possible ex- 
planation of the PAMELA positron excess, provided the boost factor in the Milky way is 
~ 3 X 10^. This interpretation is then in tension with the upper bounds for BM3 shown 
in Table [2l unless there is a mechanism which boosts the positron signal significantly 
more than the antiproton signal. 



5 Conclusions 



We have analyzed the antiproton production in the annihilation of dark matter Majorana 
particles into two leptons and one weak gauge boson. We are motivated by the possibil- 
ity of observing a gamma-ray signal from the internal electromagnetic bremsstrahlung 
process XX ~^ f fli which produces a characteristic feature in the gamma-ray spectrum 
and which may constitute a smoking-gun for dark matter detection. We have argued 
that if such a signal is observed in gamma-rays, the SU{2)l x U{1)y gauge invariance of 
the Standard Model automatically implies annihilations into two fermions and one weak 
gauge boson. Then, the fragmentation of the fermions and weak gauge bosons produces 
an antiproton flux, which is severely constrained by the PAMELA measurements of the 
cosmic antiproton-to-proton flux ratio. As a consequence, the prospects to observe a 
gamma-ray signal from internal electromagnetic bremsstrahlung at currently operating 
lACTs and at the projected CTA is constrained by the PAMELA measurements on the 
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antiproton-to-proton ratio. 

Concretely, we have studied a model where the dark matter particle is a Majorana 
fermion which couples only to the electron doublet and to an extra scalar doublet, heavier 
than the dark matter particle. We have calculated analytically and numerically the cross- 
sections for the 2—7-2 and 2—7-3 processes in all relevant scenarios, which we catalogue 
according to the level of degeneracy between the masses of the charged and the neutral 
components of the scalar doublet, and the level of degeneracy between the masses of the 
charged scalar and the dark matter particle. For all these cases, we have also calculated 
the expected antiproton-to-proton flux ratio at Earth for different dark matter halo 
profiles and different propagation models, as well as the respective upper bound on the 
cross- sect ions from requiring no excess in the PAMELA p/p data at the 95%C.L. We have 
also translated the upper bounds on the cross-section into upper bounds on the boost 
factor in the Milky Way, under the assumption that the Majorana fermion in our model 
has a thermal relic abundance which coincides with the cosmic dark matter abundance 
determined by the WMAP. 

Given the very special characteristics of this toy model, where the only source of 
antiproton production is the fragmentation of the weak gauge boson produced in the 2 — > 
3 process, our results can be interpreted as very conservative upper bounds on the cross- 
sections and boost factors in the Milky Way for more realistic models, where the dark 
matter particle also couples to quarks and weak gauge bosons. We have illustrated this 
fact calculating the upper bounds on the boost factor for some concrete benchmark points 
in the supersymmetric parameter space where the lightest neutralino constitutes the dark 
matter of the Universe. Our results indicate that in realistic models the constraints on the 
boost factors are at least one order of magnitude stronger than the constraints derived 
in this paper for the toy model. 
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Note added 

During the completion of this paper, a revised version of [l5] appeared, calculating upper 
bounds on the 2-^-3 annihilation cross-section from the PAMELA p/p data in the case 
m^o = mrj± ^ rriDM, using the corrected results from [44] . 
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Figure 13: Feynman diagrams contributing to the annihilation channel XX ~^ ^ 

A Cross-section for internal electroweak brems- 
strahlung 

The to two-to-three annihilation channel xX ~^ W~ev can proceed via radiation of a 
M/^-boson off one of the final state particles, or off the intermediate scalar particle, which 
yields three diagrams. Since x is a Majorana fermion, the two particles in the initial 
state can be exchanged, so that six diagrams contribute in total, with a relative minus 
sign for the latter ones. The differential cross-section, multiplied by the relative velocity 
is given by 



dEwdEe 16(27r)3m 



2 

DM 



-\Mi + M2 + Ms-Mi-M5-M6\^ . (A.l] 



The s-wave contribution can be obtained by setting = p^^ = (m/jjv/, 0, 0, 0). After 
expressing the neutrino momentum by the other momenta, it is easy to see that all 
remaining scalar products depend only on the energies Ew of the VT-boson and E^ of 
the e, e.g. ■ pw = fnoAtEw- In particular, from the on-shell condition pi = and 
energy-momentum conservation one obtains Pe ■ pw = ^DM{Ee + E\y) — 1m\)i^j — \M'^. 
The ly-spectrum can be obtained by integrating over the energy E^. The kinematic 
boundaries are given by i^™" = uidm — \{Ew + — M"^) and E"^"-^ = moM — 

\{Ew — ^J E"^ — M^). The result can be decomposed as 

'^da{xx W~ev)/dEw = vdawr/dEw + vdaw^/dEw ■ 
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We find the following result, in terms of x 
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(A.3) 

/i = (/i^o + /i^±)/2, 
= e^/(47r), and used 



(A.4) 



The energy of the VT-boson can vary in the range 

Mw/ruDM <x = Ew/mDM < Xmax = E'^^'^'/mDM = 1 + M^/(4m^jv^) . 

As expected, the longitudinal contribution is proportional to the square of the ri^r]^G~- 
coupling (m^± —■m'^o)/M^r to the Goldstone mode, and vanishes in the limit m^o — > m^±. 
Furthermore, the transversal contribution is proportional to Xmax~x. Thus, the transver- 
sal part of the ly-spectrum is suppressed near the endpoint, while the longitudinal part is 
not. This behaviour is expected from angular momentum conservation, as was explained 
in Section [21 Note that the charge-conjugated process XX ~^ W^ev produces an identical 
spectrum of in addition. 

The spectrum of Z-bosons produced in the electroweak IB annihilation processes 
XX ~^ and xX ~^ can be obtained from the above formulae by replacing 

nirp m^± or m^± i— )■ m^o, respectively, as well as My/ i-> Mz- Furthermore, one has to 
multiply by 2s^ tan^(26'i4/) for xX ~^ ^ee and by l/(2c^) for xX ~^ Zuu. 

Finally, we note that one can obtain the photon spectrum produced in the electro- 
magnetic IB process xX ~^ 7^6 by replacing m^o H- m^±, then taking the limit My/ — )■ 0, 
and multiplying by 2s^. The resulting formula coincides with Eq. 9 of Ref. [8] (including 
the prefactor). We have checked that in the limit m^o — > m^± our result for the W spec- 
trum is also consistent with the result reported in Ref. [43] (also including prefactors), 
and with the result of Ref. up to an overall factor of two. 
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The total cross-sections can be obtained by integrating over the energy of the gauge 
boson. Here, we just present the result obtained in the limit m^i ^ ttidm, when keeping 
m^± — rriio fixed (for our analysis, we use the full expressions). 



^4 



+ 5A/i2(l + lOr + r^)] - 30r(l + r)(2r - A/i^) ln(r)| , (A.5) 
avixx^W^eu) = "^'"C , Z^"":" ~ ^'^^^ |l - 8r + Sr^ - - 12rMn(r)| , 

where r = M^/(4m^^j). The cross-sections for XX ~^ Zee and XX ~^ Zuu can be 
obtained by the same manipulations as described above for the spectra. In the limit 
m^i ^ rrioAi, the total cross-section for electromagnetic IB is, for comparison, given by 
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